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Polarisation of high-energy emission in a pulsar striped wind 
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Abstract. Recent observations of the polarisation of the 
optical pulses from the Crab pulsar motivated detailed 
comparative studies of the emission predicted by the polar 
cap, the outer gap and the two-pole caustics models. 

In this work, we study the polarisation properties 
of the synchrotron emission emanating from the striped 
wind model. We use an explicit asymptotic solution for 
the large-scale field structure related to the oblique split 
monopole and valid for the case of an ultra-relativistic 
plasma. This is combined with a crude model for the emis- 
sivity of the striped wind and of the magnetic field within 
the dissipating stripes themselves. We calculate the polar- 
isation properties of the high-energy pulsed emission and 
compare our results with optical observations of the Crab 
pulsar. The resulting radiation is linearly polarised. In the 
off-pulse region, the electric vector lies in the direction of 
the projection on the sky of the rotation axis of the pulsar, 
in good agreement with the data. Other properties such 
as a reduced degree of polarisation and a characteristic 
sweep of the polarisation angle within the pulses arc also 
reproduced. 



1. Introduction 

The high-energy, pulsed emission from rotating magne- 
tised neutron stars is usually explained in the framework 
of either the polar cap or the outer gap models. Although 
the existence of such gaps is plausible (Petri et al. 2002), 
these models still suffer from the lack of a self-consistent 
solution for the pulsar magnetosphere. Nevertheless, re- 
cent observations of the polarisation of the optical pulses 



from the Crab (Kanbach et al. 2003) motivated detailed 
comparative studies of the emission predicted by the po- 
lar cap, the outer gap and the two-pole caustics mod- 



els (Dyks et al. 2004) In all of these models, the radia- 
tion is produced within the light cylinder. However, the 
pulse profile is determined by the assumed geometry of 
the magnetic field and the location of the gaps. None of 
these models is able to fit the optical polarisation proper- 
ties of the Crab pulsar. 



An alternative site for the production of pulsed radi- 

based on 



ation has been investigated (Kirk et al. 2002) 



the idea of a striped pulsar wind, originally introduced 



by Coroniti (1990) and Michel (1994) Emission from the 



striped wind originates outside the light cylinder and rel- 
ativistic beaming effects are responsible for the phase co- 
herence of the synchrotron radiation. A strength of this 
model is that the geometry of the magnetic field, which is 
the key property determining the polarisation properties 
of the emission, is relatively well-known. 

In this work, we use an explicit asymptotic solution for 
the large-scale field structure related to the oblique split 
monopole and valid for the case of an ultra-relativistic 



plasma (Bogovalov 1999) This is combined with a crude 



model for the emissivity of the striped wind and of the 
magnetic field within the dissipating stripes themselves. 
We calculate the polarisation properties of the high-energy 
pulsed emission and compare our results with optical ob- 
servations of the Crab pulsar. 

2. Stokes parameters 

Our magnetic field model is based on the asymptotic so- 
lution of the split monopole for the oblique rotator, valid 
for r>rL, and modified to take account of a finite width 
of the current sheet. We add a small meridional compo- 
nent in this sheet in order to prevent the magnetic field 
from becoming identically zero. The 3-dimensional geom- 
etry of the current sheet is shown in figure [TJ In spherical 
polar coordinates (r, 9, ip) centered on the star and with 
axis along the rotation axis, the radial field is small and 
neglected, B r ~ B^r^/r 2 , while the other components 
are: 



B- 
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Fig. 1. 3D structure of the current sheet in the striped 
wind. The rotating neutron star is located at the origin of 
the coordinate system. 

Here, _Bl is a fiducial magnetic field strength, v is the 
(radial) speed of the wind, £7* = c/r-^ is the angular 
velocity of the pulsar, with tl the radius of the light 
cylinder, a is the angle between the magnetic and rota- 
tion axes, &i 2 are parameters controlling the magnitude 
of the meridional field in the two current sheets present 
in one wavelength, and Ag i¥ , are parameters quantifying 
the sheet thickness. The functional form of B v is moti- 
vated by exact equilibria of the planar rclativistic current 



sheet (Kirk & Skjaeraasen 2003) However, in these equi 



libria the Bg component, which has an important influence 
on the polarisation sweeps, is arbitrary. The Bg we adopt 
corresponds to a small circularly polarised component of 
the pulsar wind wave, such as is expected if the sheets are 
formed by the migration of particles within the wave, as 



described qualitatively by Michel (1971) 



For the particle distribution, wc adopt an isotropic 
electron/positron distribution given by 



N(E,p,r,t) = K(r,t)E~ 



(5) 



where K(r, t) is related to the number density of emitting 
particles. The radial motion of the wind imposes an overall 
1/r 2 dependence on this quantity, which is further mod- 
ulated because the energization occurs primarily in the 
current sheet. The precise value in each sheet is chosen to 
fit the observed intensity of each sub-pulse. In addition, a 
small dc component is added, giving the off-pulse inten- 
sity. For the emissivity, we use the standard expressions 
for incoherent synchrotron radiation of ultra-relativistic 
particles. We assume the emission commences when the 
wind crosses the surface r = ro 3> t"l. 

The calculation of the Stokes parameters as measured 
in the observer frame involves simply integrating the emis- 
sivity over the wind. For an observer, at time t bsj they 



are given by the following integrals: 



Quj ^ (tobs) = 

u.. I 



■OOf7T />2 77 

So(r, tret, 



P+7/3 
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cos (2 x) 
sin(2x) 



d 3 r(6) 



where the retarded time is given by t rct = i b s + n ■ r/c 
and n is a unit vector along the line of sight from the 
pulsar to the observer. In this approximation the circular 
polarisation vanishes: V = 0. The function sq is defined 
by: 



s (r,t) = nK{r,t) 



V 2 
U! 2 



|v/l-(2?n.6)2 



(7) 



where u> is the angular frequency of the emitted radiation, 
and At is a constant factor that depends only on the nature 
of the radiating particles (charge q and mass m) and the 
power law index p of their distribution: 
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with Teu the Euler gamma function and V the Dopplcr 
boosting factor V = 1/r (1 — (3 ■ n). The direction of the 
local magnetic field in the observer's frame is given by the 
unit vector b and the simplifying assumption has been 
made that this field has no component in the direction of 
the plasma velocity: b ■ (3 = 0. (In this case the magnetic 
field transformation from the rest frame B' to the observer 
frame B is just B' = B/T and, thus, its direction remains 
unchanged.) The angle x measures the inclination of the 
local electric field with respect to the projection of the 
pulsar's rotation axis on the plane of the sky as seen in 
the observer's frame. The degree of linear polarisation is 
defined by 



n 



y/Q 2 + U 2 



(9) 



The corresponding polarisation angle, defined as the posi- 
tion angle between the electric field vector at the observer 
and the projection of the pulsar's rotation axis on the 
plane of the sky is 



1 fU 
X = — arctan | — 



(10) 



3. Results 



The upper left panel of Fig. [5] shows the intensity (Stokes 
parameter /) computed using our smoothed profile with 
Ag = 1, A^, = 5, b\ = 0.1 and 62 = 0.08 for each sub- 
pulse. The electron density is 



K 



r 2 (1- 0.6 rjg) 
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Fig. 2. Light curve of intensity, degree of polarisation and 
position angle of the pulsed synchrotron emission obtained 
for our model and measurements of these quantities for 
the Crab pulsar. Models with Lorentz factor T = 20 (solid 
red) and 50 (dotted cyan) are shown. The bottom panels 
show the dependence on phase of the assumed magnetic 
field components and the particle density in the comoving 
frame. 



where the parameter e = 0.05 sets the minimum electron 
density between the current sheets (in normalized units). 
The denominator (1 — 0.6 r]g) introduces an asymmetry 
in the relative pulse peak intensity. The variation of the 
magnetic field and the particle density along the line of 
sight, are shown in the bottom panels of Fig. [21 

The upper panels of Fig. [5] show the results of our 
computations on the left and the corresponding observed 



quantities (Kanbach et al. 2003) on the right. Comparison 
with the upper right-hand panel shows that the model 
reproduces the observed pulse profile quite accurately. 
In this example, we adopt an obliquity a = 60° and 
an inclination of the rotation axis to the line of sight 
£ = arccos(n ■ f2*/|fi*|) = 60° and set the radius at 
which emission switches on to be r$ = 30rL- The posi- 
tion angle of the projection of the pulsar's rotation axis 
is set to 124° |(Ng fc Romani 2004)| The electron power 



law index p = 2, as suggested by the relatively flat spec- 
trum displayed by the pulsed emission between optical 



and gamma-ray frequencies (Kanbach 1998) Results are 
shown for two values of the Lorentz factor of the wind: 
L = 20 (solid line) and P = 50 (dotted line). For con- 
venience, the maximum intensity is normalized to unity. 
The timcscalc is expressed in terms of pulse phase, corre- 
sponding to the initial time t = and 1 to a full revolution 
of the neutron star and thus to one period = 2n/£l*. 



The degree of polarisation is shown in the two middle 
panels of Fig. [21 According to our computations (left-hand 
panel) this displays a steady rise in the initial off-pulse 
phase, that steepens rapidly as the first pulse arrives. Dur- 
ing the pulse phase itself, the polarisation shrinks down 
to about 10%. Theoretically the maximum possible degree 
of polarisation is closely related to the index p of the par- 
ticle spectrum. In the most favorable case of a uniform 
magnetic field, it is given by 



p + i 

p + 7/3 



(12) 



However, in the curved magnetic field lines of the wind, 
contributions of electrons from different regions have dif- 
ferent polarisation angles. Consequently, they depolarise 
the overall result when superposed. We therefore expect a 
degree of polarisation that is at most II max . For the exam- 
ple shown in figure [21 n max (p = 2) = 69.2%, well above 
the computed value, which peaks at 52%. 

The lower panels of Fig. [21 show the polarisation an- 
gle, measured against celestial North and increasing from 
North to East. Our model predicts this angle relative to 
the projection of the rotation axis of the neutron star on 
the sky, which we take to lie at a position angle of 124°, 



following the analysis of (Ng & Romani 2004) In the off- 
pulse stage, the electric vector of our model predictions 
lies almost exactly in this direction, since it is fixed by the 
orientation of the dominant toroidal component B v of the 
magnetic field. In the rising phase of the first pulse, B v 
decreases, whereas Bg increases, causing the polarisation 
angle to rotate from its off pulse value by about 50°, for 
the chosen parameters. However, for a weaker Bg contri- 
bution, as in the second pulse, the swing decreases. This 
effect can also be caused by a relatively large beaming an- 
gle, (i.e., low Lorentz factor wind). The basic reason is that 
contributions from particles well away from the sheet cen- 
ter are then mixed into the pulse, partially canceling the 
contribution of the particles in the center of the current 
sheet, which favor x = 124° ±90°, and enforcing x — 124°. 
On the other hand, a very high value of the Lorentz factor 
or large values of reduce the off-center contribution, 
leading, ultimately, to the maximum possible 90° sweep 
between off-pulse (-B v -dominated) and center-pulse (Bg- 
dominated) polarisations, followed by another 90° sweep 
in the same sense when returning to the off-pulse. Thus, 
in general, in the middle of each pulse, the polarisation 
angle is either nearly parallel to the projection of the ro- 
tation axis, or nearly perpendicular to it, depending on 
the strength of Bg and on T. This interpretation is con- 
firmed by computations with T = 50 that show a larger 
sweep, as shown in Fig. [2] 

The optical polarisation measurements suggest that in 
the centre of the pulses the position angle is close to 124°. 
In the declining phase of both pulses, the angle reaches 
a maximum before returning to the off-pulse orientation. 
Note that in both cases the swing starts in the same di- 
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rection, (counterclockwise in figure [2|). This is determined 
by the rotational behavior of the Bg component, imply- 
ing that this changes sign between adjacent sheets, as in 
Eq. (|4]). The observed off-pulse position angle is closely 
aligned with the projection of the rotation axis of the pul- 
sar, in accordance with the model predictions. 

In addition to models aimed at providing a framework 
for the interpretation of the emission of the Crab pul- 
sar, we have performed several calculations with different 
Lorentz factors T, injection spectrum of relativistic elec- 
trons p and inclinations of the line of sight £. The general 
characteristics of the results are: For low Lorentz factors, 
independent of p and £, the relativistic beaming becomes 
weaker and the pulsed emission is less pronounced, be- 
cause the observer receives radiation from almost the en- 
tire wind. For instance, taking T = 2 and p = 2 or 3, the 
average degree of polarisation does not exceed 20 % and 
the swing in the polarisation angle is less than 30°. For 
high Lorentz factors V > 50, the strong beaming effect 
means that the observer sees only a small conical frac- 
tion of the wind. The width of the pulses is then closely 
related to the thickness of the transition layer. The de- 
gree of linear polarisation flattens in the off-pulse emis- 
sion while it shows a sharp increase followed by a steep 
decrease during the pulses. Due to the very strong beam- 
ing effect, only a tiny part of the wind directed along the 
line of sight will radiate towards the observer. In the off- 
pulse phase, the polarisation angle is then dictated solely 
by the B v component "attached" to the line of sight, and 
the degree of linear polarisation remains almost constant 
in time. For very high Lorentz factors, the behavior of 
polarisation angle and degree remain similar to those of 
Fig. [2j with perfect alignment between polarisation direc- 
tion (electric vector) and the projection of the pulsar's 
rotation axis on the plane of the sky in the off-pulse phase 
and two consecutive polarisation angle sweeps of 90° in 
the same sense during the off-pulse to center-pulse and 
center-pulse to off-pulse transitions. This mirrors the fact 
that emission comes only from a narrow cone about the 
line of sight of half opening angle 6 « 1/T. 

For given values of T and £, the particle spectral index 
p affects only the average degree of polarisation degree but 
not the light curve nor the polarisation angle. For example, 
taking r = 10 and £ = 60°, a spectral index of p = 2 leads 
to an average polarisation of ft = 19.2% whereas for p = 3 
it leads to ft = 30.8%. 



them on the basis of observations. However, the geometry 
of the magnetic field, which is the crucial factor deter- 
mining the polarisation properties, is constrained in the 
striped model to be close to that of the analytic asymp- 
totic solution of the split monopole. We have therefore 
presented detailed computations of the polarisation prop- 
erties of the pulses expected in this scenario. These pos- 
sess the characteristic property, unique amongst currently 
discussed models, that the electric vector of the off-pulse 
emission is aligned with the projection of the pulsar's ro- 
tation axis on the plane of the sky. This is in striking 
agreement with recent observations of the Crab pulsar. In 
addition the striped wind scenario naturally incorporates 
features of the phase-dependent properties of the polari- 
sation angle, degree of polarisation and intensity that are 
also seen in the data. This underlines the need to develop 
the model further, in order to confront high-energy obser- 
vations of the Crab and other pulsars. In particular, the 
manner in which magnetic energy is released into parti- 
cles in the current sheet remains poorly understood and 
the link between the asymptotic magnetic field structure 
and the pulsar magnetosphere is obscure. 
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4. Conclusions 

In the striped wind model, the high energy (infra-red to 
gamma-ray) emission of pulsars arises from outside the 
light cylinder. It provides an alternative to the more in- 
tensively studied gap models. They all contain essentially 
arbitrary assumptions concerning the configuration of the 
emission region and the distribution function of the emit- 
ting particles, rendering it difficult to distinguish between 



